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1, INTRODUCTION 

The world electrical power demand has increased about 22.9% from the year 2005 to the year 2016 
due to rapid growth in industrial and commercial activities [1]. Due to this situation, several options have 
been considered to meet the future energy demand for the existing power system. One of the options is the 
integration of renewable energy such as the solar energy into the existing grid. Distributed Generation 
Photovoltaic (DGPV) is the preferable source commonly implemented in power system. The role of DGPV is 
mainly to provide the active power and no reactive power generated to the system. This solution not only 
able to meet the increasing power demand but also can further improve the power losses and the system 
stability [2]. Effect of DGPV on the distribution [3],[4] and transmission systems [5] have been an interesting 
subject for many researchers. Most of the research efforts are conducted for location and sizing of DGPV to 
satisfy the technical benefits such as loss minimization, voltage stability enhancement and maximum 
loadability increment. Among these, some of the researchers focused only to minimize the power losses as 
the objective function [6]-[8]. In other studies, the optimal location and sizing of DG were studied to focus 
only on the voltage stability improvement as the objective function [9],[10]. 

In the last few years, various techniques have been developed to find the optimal location and size 
of DG. These techniques can be categorized as the analytical methods and meta-heuristic methods. Several 
analytical approaches are proposed for DG allocation to minimize the power losses [11],[12]. The review of 
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most analytical methods for DG allocation is discussed in [13]. For the same purpose, evolutionary algorithm 
(EA) techniques have also been applied for single or multi DGs location and sizing for various objectives. 
EA techniques are preferable due to its computational times compared to analytical techniques. EA 
techniques start with random initialization of the population followed by the evolvement of the population 
across several generations. In each generation, fit individuals are selected to become parent individuals. Then 
these individuals undergo the mutation process to produce offspring individuals. 

In evolutionary algorithm, there are two important criteria which need to be considered in order to 
enhance the performance of the algorithm namely the exploration and the exploitation phase [14],[15]. 
Exploration space refers to the ability of the algorithm to search for a solution in the whole region of search 
space. Meanwhile, exploration specifies the convergence towards the best optimal solution in the exploration 
space. There are many techniques used to improve these phases for example gradient descent, random walk 
and local search. However, currently, many researchers have considered using chaos theory as one of the 
approaches to improve the performance of evolutionary algorithm [16]-[18]. From the previous studies, many 
metaheuristic algorithms have been integrated with a chaotic map in order to improve the performance. In 
[18], the chaotic local search is utilized in Differential Evolution (DE) and has been tested on 13 classical test 
functions. The result showed significant improvement in exploitation phase as compared to the traditional 
DE. Peng Lu et al [19] also concluded that chaotic behaviour is able to enhance the performance of 
differential bee colony optimization to solve the economic dispatch problem. The results of these studies 
provide the proof of how successful chaos theory is in improving the evolutionary algorithm. In this study, 
the chaotic mapping is used to enhance the performance of Chaotic Mutation Immune Evolutionary 
Algorithm (CMIEP) as the local search technique for optimal DGPV allocation. Comparative studies were 
performed with respect to Evolutionary Programming (EP). Results had indicated that CMIEP with chaotic 
local search outperformed EP in terms of accuracy of FVS/ and transmission losses. 


2. PROBLEM FORMULATION 
In this study, two single objective functions are considered and optimized separately while 
satisfying system equality and inequality constraints. 


2.1. Objective functions 
The objective of the optimal location and sizing of DGPV is to minimize two objective functions 
namely the active power loss and the voltage stability index separately. 


2.2. Minimization of Voltage Stability Index 

A line based voltage stability index, FVS7 developed by I. Musirin et al. [20] is used to measure the 
closeness of the system to voltage collapse. The FVS/ formulation was derived from a voltage quadratic 
equation on a two-bus system and defined by the following equation: 





(1) 


is line impedance 

is line reactance 

is reactive power at the receiving end 
is voltage at the sending end 

A(X) is the first objective function 


aie re} 


2.3. Minimization Active Power Loss 
The total real power loss, Piss in the transmission lines can be expressed as in equation (2): 


fo (X) = P55; forVvienr 
F,(X)=min(f, (x) ” 


nr is the number of transmission lines. 
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2.4. Constraints 
The objective functions are subjected to the following constraints: 
1) The generating capacity 


PGi S Poa < DG i VIE N (3) 


where PG, and PZ, are the minimum and the maximum output of DGPV respectively and 7 is the total 
bus number. The bus voltage constraint is defined as follows: 


Vmin < Vi < Vmax V1 = N (4) 


where Vin and Vinax are the lower and the upper bound of bus voltage limit respectively and V; is the voltage 
magnitude at bus i for all the NV bus. The power balance constraint is shown in equation (5): 


(PGi-Ppit Poci)= Ploss ViEN (5) 


where Pg;i,Ppi,. PpGjand Pioss are the active power of bus generator, active load and active power 
losses respectively. 7 1s the total bus number. 


3. CMIEP FOR DGPV PLACEMENT 

For a transmission network, load flow analysis is carried out and FVS/ or loss value is computed for 
each line using Equation (1) and Equation (2) respectively. The CMIEP algorithm is used for finding the 
optimum size of DGPV at an optimum location based on a minimum total power loss, with constraints given 
in Equation (3) to (5). In this study, the 30 Bus IEEE RTS is used as the test system. The complete flow chart 
for DGPV allocation and sizing 1s represented in Figure 1. 


3.1. Chaotic Local Search 

To improve the search capability and for achieving global optima solution of DGPV location and 
sizing, chaotic dynamics is incorporated into CMIEP. A chaotic function known as Piecewise Linear is 
employed and the equation is defined as in Equation (6) [21],[22]: 


c,/ p c, €( 

a (c; - p)/(0.5- p) c, €( 
ae (1- p-c,)/(0.5- p) c, €( 
l-c,/p c, €( 


(6) 


Cc, is the chaotic variable that is influenced by the value of control parameter, p. The PWLCM exhibits 


t 
chaotic dynamics in (0,1) when control parameter, p € (0,0.5) U (0. 5,1). 
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Figure 1. Flowchart for DGPV Allocation Using CMIEP 


Step 1: Setting the iteration, t=0, global best decisions, Xyoor s(t) i—1,2...,V_ ,fitness value from the 


optimization F (x! ry ) initial value of search space for each of the variable, 7; as in Equation (7): 


1()-(2ZB) a-12.1 


(7) 
and chaotic variables, c,(t)= rand(1, x) where UB and LB are the upper and lower boundary of 


searching space for decision variables x. and V is the total decision variables to be solve respectively. 
Step 2: Determine the chaotic variables ¢ (+1) for the next iteration using the chaotic equation in 


equation (6) according toc, (t). 


Step 3: Determine the decision variable U : by converting the chaotic variables C, (t + 1) using 


the following equation : 


l 


Uj = 4, (+1) =x, (t) +(e, (¢4+1)-05)* 7, (2) (8) 


Step 4: Calculate the new fitness Flu! ) 


Step 5: Evaluating the fitness value of the new solution, Flu‘) with the optimize fitness value 


F (x - r ) using the evaluation step as Figure 2: 
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/* Evaluation step * / 

t=0 

while t <50 

if F(U; )SF (%pen;) then 
Gea. 
xh xe 


else 


t t 
Dias = Nest, 
end if 
fats] 


end while 





Figure 2. Evaluation Algorithm 


Step 6: Update the search radius for each decision variables, x; as in Equation (9): 


l, (1 + 1) =F, (1) * rand (0,1) (9) 


Step 7: If maximum iteration is reached, display the output of chaotic local search. Otherwise, go 
back to Step 2. 


4. SIMULATION RESULTS 

The chaotic variable is initialized by the rand function whose advantage is already addressed in the 
previous section. IEEE 30 bus RTS is used as the test system to verify the feasibility and robustness of the 
proposed CMIEP with chaotic local search. Two DGPV installations are used in this study to find optimal 
location and size of DGPV. There are three cases considered in this paper to monitor the capability of the 
optimization technique. The simulation result obtained using CMIEP is then compared with the result 
obtained by EP to prove the effectiveness of the proposed algorithm. 


4.1. Base Case 

The best results obtained by the implementation of the two objectives separately in base case 
condition is shown in Table 1. The results of post-installation of DGPV have been compared with the results 
from pre-DGPYV installation. The FVS/ of post-installation using CMIEP has been reduced to 32 % from the 
pre-installation FVS7. Meanwhile, by using EP the FVSI has been reduced to 25% of the pre-installation 
value. The loss shows the reduction of about 67% from the pre-installation value by using CMIEP and about 
65% by using EP. Therefore, it can be seen that the proposed algorithm is capable of finding better solutions 
for each objective as compared to EP. The convergence characteristic for this case is presented in Figure 3. 
Results indicated that CMIEP has better fitness and convergence rate compared to EP. 


Table 1. Comparison of Objective Functions for 2 DGPV Installation 


ae Post-Installation 
Objective Function Pre-installation 


CMIEP EP 
FVSI 0.2037 0.1381 0.1519 
Loss (MW) 17.58 5.68 5.99 
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Figure 3. (a) FVS/ and (b) Loss Convergence Characteristics of CMIEP and EP with Two Units DGPV 


4.2. FVSI Minimization 

Results for FVS7 minimization using CMIEP when load bus 29 is subjected to load variation 1s 
tabulated in Table 2. The increment of loading condition showed a significant increase on the maximum 
FVSI of the system. However, after the DGPV installation, FVS/7 has been reduced to28%, 6% and 9% for 


load variation of 10, 20 and 30 MVAR respectively. The location and sizing of DGPV to achieve the 
improvement of F'VS/ can be referred to the same table. 


Table 2. FVST Minimization when Load Variation was Subjected to Bus 29 


Loading Post-Installation using CMIEP 
Condition Pre-Installation 
DGPV DGPV % FVSI 
Qa2o(MV AR) FVSI Location Size (MW) FVSI Reduction 
10 0.2111 17 29 42.89 42.90 0.1503 28.78 
20 0.3573 30 24 13.18 56.29 0.3359 5.99 
30 0.5987 30 24 15.96 56.61 0.5449 8.99 


4.3. Transmission Loss Minimization 


Table 3 tabulates the result for DGPV optimal location and sizing using CMIEP. With the same 
loading conditions as the previous case, DGPV can also reduce the transmission losses of the system. For 
instance, at loading condition of 20 MVAR the transmission losses reduced to 7.19 MW which 62.91% 
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reduction from the pre-installation loss. The optimal location and sizing optimized by CMIEP technique are 
buses 21 and 7 with 59.68 MW and 61.66 MW respectively. From the table, it is clearly showed that bus 12 
and bus 7 are the optimal location for all the loading condition. 


Table 3. Transmission Loss Minimization when Load Variation was Subjected to Bus 29 


Lesdine Condition Pre-Installation Post-Installation using CMIEP 
Loss DGPV DGPV Loss % Loss 
Quo (MV AR) (MW) Location Size (MW) (MW) reduction 
10 18.12 21 7 59.68 64.01 6.00 66.87 
20 19.39 21 fi 56.37 61.66 7.19 62.91 
30 22.44 21 d 60.65 50.56 10.10 55.00 


The results of minimization of FVS/ and minimization of losses for CMIEP and EP are compared 
with the pre-installation value and shown in Figure 4. In Figure 4(a), both CMIEP and EP are comparable at 
load variation of 10 MVAR and 20 MVAR. However, at load variation of 30 MVAR, CMIEP outperformed 
EP with 4% difference of post-installation loss reduction. In Figure 4(b), it is clear to mention that CMIEP 
outperformed EP in all cases to determine the optimal location and size of DGPV for FVS/ reduction. 
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20.00 


15.00 


oss (MW) 


. 


I 


10.00 


5.00 





10 20 30 
Load at Bus % Qd (MVAR) Load at Bus 29, Qd (MVAR) 


(a) (b) 


0.00 


Figure 4: (a) Loss (b) FVS/ with Load Variation Subjected to Bus 29 


5. CONCLUSION 

This study proposed a modified version of a pre-developed Chaotic Mutation Immune Evolutionary 
Programming (CMIEP). Chaotic local search has been added into the original CMIEP algorithm. The 
inclusion of chaotic local search managed to achieve a better optimal solution for location and sizing of 
DGPYV in the transmission system. In the proposed placement scheme, the transmission loss and FVSI/ value 
are treated as the fitness equations and optimized separately. From the study, the optimal value of DGPV 
units for different loading levels are changed as load changes. The results also revealed that the utilization of 
DGPV into the transmission system reduces the total of active power losses and effectively improve the 
voltage stability of the system. A comparative study, also showed that the proposed CMIEP algorithm 
outperformed EP exhibited by a better reduction in FVS/ values and lower loss values, for all cases. In the 
future, CMIEP can be used to perform multi-objective functions taking FVS/ and losses as the objective 
functions. 
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